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Abstract 

We present a one-parameter family of periodic packings of regular 
tetrahedra, with the packing fraction 100/117 ~ 0.8547, that are simple 
in the sense that they are transitive and their repeating units involve only 
four tetrahedra. The construction of the packings was inspired from re- 
sults of a numerical search that yielded a similar packing. We present an 
analytic construction of the packings and a description of their proper- 
ties. We also present a transitive packing with a repeating unit of two 
tetrahedra and a packing fraction 139 \$g V ^ ~ 0.7194. 

1 Introduction 

The optimization problem of packing regular tetrahedra densely in space has 
seen invigorated interest over the last few years [3j HJ [6l [7] • This interest has 
helped drive up the packing fraction of the densest-known such packings from 
0.7174 in 2006 [3J to 0.8503 [7] most recently (see Table 3). These improved 
packing fractions have been obtained from more and more complex packings, 
with larger and larger repeating units. This trend has led some to conjecture 
that the densest packing of tetrahedra might have inherent disorder [BJ. The 
more restrictive problem of packing tetrahedra transitively — that is, so that 
all tetrahedra in the packing are equivalent (a more rigorous definition is given 
below) — has been less extensively studied and the densest previously-reported 
transitive packing of regular tetrahedra fills only 2/3 of space |3J. Here we 
present a one-parameter family of transitive but dense packings of tetrahedra 
with the packing fraction 100/117 ks 0.8547. 

The discovery of this family of dense packings was inspired by the results 
of a numerical search, which yielded a dense packing with similar structural 
properties to the packing we present. The numerical method used was adapted 
from the divide and concur approach to constraint satisfaction problems [2\ . The 
divide and concur formalism enables us to set up an efficient search through the 
parameter space consisting of the positions and orientations of tetrahedra inside 
the repeating unit and the translation vectors governing its lattice repetition, 
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fundamental 
tetrahedron 


T = convjr, 1 i = 1,2,3,4} 
r =Ia-lb + §c 

1 28 c? 39 

r 2 = ia - 

r 3 = T4 a +To b +n c 
r 4 = fa+ib-lc 


other tetrahedra 
in the unit cell 


Ti = R 2 (T ) = conv{|(r 2 + r 3 + r 4 ) - ri,r 2 ,r 3 ,r 4 } 
T 2 = /(T ) = -T 
T 3 = /(TO = -Ti 


space group 
generators 


translations by a, b, = |b + + ia 
I = inversion about 

i?2 = two- fold rotation about {|a + th \ t G K} 


packing fraction 


100/117 


coordinate basis for which 
tetrahedra are regular 
with unit edge length 


a = (2^7/5,0,0) 

b= (0,^3/2,0) 

c = (0,0,13^/3/14/5) 



Table 1: The construction of the dimer-double- lattice family of packings in 
terms of the parameter 29/56 < x < 9/14 in a general monoclinic coordinate 
basis (a • b = b ■ c = 0) . The packing is generated starting from the fundamental 
tetrahedron by the action of the space group. A packing of regular tetrahedra is 
obtained when the general monoclinic coordinate basis reduces to the specified 
orthogonal coordinate basis. 

subject to the constraint that no two tetrahedra overlap. The dynamics involved 
in the divide and concur search are highly non-physical, which might explain 
why our method was able to discover this dense packing, while earlier methods 
involving more physical dynamics were not [SJ El [Z] • I n this note we present only 
the analytically constructed packing without a full explication of the numerical 
method, which will be forthcoming. 

2 One-parameter family of dimer-double-lattice 
packings 

The first set of packings we report on are naturally described as double lattices of 
bipyramidal dimers. A double lattice is the union of two Bravais lattices related 
to each other by an inversion operation about some point. In [T], Kuperberg and 
Kuperberg used the idea of a double lattice in the Euclidean plane to show that 
any planar convex body can be packed in an arrangement with a packing fraction 
no smaller than y/3/2. We naturally extend the idea of the double lattice to the 
three-dimensional Euclidean space. The repeating unit of one constituent lattice 
is a bipyramidal dimer: two regular tetrahedra sharing a common face. Two 
of these dimers with mutually-inverted orientation — a Kuper-pair — form the 
repeating unit of the double lattice, which thus has four tetrahedra of distinct 
orientations. We state the existence of the packings as Theorem 1. 
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Theorem 1. There exists a one-parameter family of packings of regular tetra- 
hedra, each having packing fraction 110/117. These packings are periodic, with 
each unit cell of the lattice containing two bipyramidal dimers. The group of 
isometries leaving each packing invariant is a crystallographic space group of 
type C2/c (following the classification and notation of \1 and acts transi- 
tively on the individual tetrahedra of the packing. 

Proof. We construct each packing by acting on a single regular tetrahedron with 
a group of isometries. For a coordinate basis we use three pair-wise orthogonal 
vectors a, b, and c, of norms |a| = 2a/7/5, |b| = y/3/2, and |c| = 13y/3/U/5. 
The initial tetrahedron To = conv{rj | i = 1,2,3,4} is the convex hull of four 
vertices whose coordinates are given in Table 1; it is a regular tetrahedron of unit 
edge length. The group of isometries is the group of crystallographic type C2/c 
generated by the translations a, b, and d x = ^b+ ^c + xa (29/56 < x < 9/14), 
by the inversion about the point 0, and by the rotation by 180 degrees about the 
axis {ja + tb \teR} [11]. This space group has a point group of order 4 and 
its translations generate a centered monoclinic point lattice. The construction 
is summarized in Table 1. 

By construction, each tetrahedron in the packing is the image of To under 
an isometry in the group. We have immediately then that all tetrahedra are 
congruent with To, that the packing is invariant under the action of the group, 
and that the group acts transitively on individual tetrahedra. As the tetrahedra 
divide into four orbits of the lattice translations, the packing fraction is given 

4vol(T ) 4(a/2/12) _ 100 

~ |det([a,b,d a ])| ~ |a||b||c|/2 ~ 117 

All that is left then to prove the theorem is to verify that the arrangement 
of tetrahedra thus constructed is indeed a packing. As the packing is transi- 
tive (in the sense that its symmetry group acts transitively on the constituent 
tetrahedra), it is enough to check that one tetrahedron, To, does not overlap 
any other tetrahedron. By means of a closest-lattice-point algorithm such as 
the one in [TO], we generate a list of all tetrahedra whose centroid is, for any 
29/56 < x < 9/14, at a distance less than y/3/2 from the centroid of T . There 
are 46 such tetrahedra. All other tetrahedra have circumspheres which do not 
intersect the circumsphere of To, and therefore they do not intersect To. For 
each tetrahedron in the list we can establish the existence of a separating plane 
separating it from To. Two tetrahedra have no overlap if and only if they are 
separated by a plane, and moreover, such a plane always exists which passes 
through three of the eight vertices of the two tetrahedra. By exhaustively veri- 
fying that one of the finitely many planes that pass through three of the vertices 
separates the two tetrahedra, we establish that each tetrahedron in the list can 
be separated from Tq. □ 

Note that the above construction, which uses a specific orthogonal coordinate 
basis {a, b, c}, is a special case of a general family of packings of non- regular 
tetrahedra that can be obtained using the same construction, but with a general 
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Figure 1: Small portions of one layer and three stacked layers in the dimer- 
double- lattice packing given by x = 4/7. 

monoclinic coordinate basis (a ■ b = b ■ c = 0). Each of these more general 
packings is an image of a packing in the original family under an affine map 
from a three-parameter family (not counting pure dilation). As the lack of 
overlap between tetrahedra and the packing fraction are both affine-invariant, 
these are also transitive packings of the same packing fraction. 

By the construction of the double lattice, there is an inversion center that 
sends one lattice of dimers into the other. Note that a lattice translation com- 
posed with an inversion about a point corresponds to an inversion about a point 
related to the original inversion center by half the lattice vector. It follows then 
that in any primitive unit cell of the lattice, there are eight such inversion 
centers. These eight inversion centers form the vertices of a parallelepiped one- 
eighth the volume of the primitive unit cell of the lattice. This parallelepiped 
is the equivalent of the "extensive parallelogram" described in [1] whose vertices 
are the inversion points that generate the double lattice. As in [Tj, the paral- 
lelepiped is inscribed in the body being packed — the bipyramidal dimer in our 
case. 

An interesting feature of the packing is the presence of the free parameter 
x. The effect of a change in x is to slide fixed layers of the packing relative 
one another along the a-direction. These layers are the layer generated by the 
translations a and b from the four tetrahedra of the primitive unit cell and the 
layers parallel to it. The construction yields a valid packing when the small 
protrusions in one layer are staggered to fit into small gaps in the neighboring 
layer, which is the case for all 29/56 < x < 9/14 (mod 1). Within this range, 
each layer can slide against the neighboring layer without changing the spacing 
between the two or creating collisions. It is possible then to obtain equally dense, 
non-transitive packings by staggering consecutive layers arbitrarily within the 
allowed range. 

We describe next the contacts formed by each dimer in the packing, and they 
are illustrated in Figures 1 and 2. Each of the eight vertices of the inscribed par- 
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Figure 2: The contacts on the surface of a dimer shown on a net diagram for 
x = 29/56 (left) and x = 9/14 (right): the face- face contacts (gray), whose 
centers (black dots) lie on inversion centers, four of which are fixed and four 
of which move as a function of x; the four point contacts made regardless of 
the value of x (blue squares), all lying on two- fold axes; the four point contacts 
formed only for x = 29/56 (purple dots); and the four point contacts formed 
only for x — 9/14 (green dots). 



allelepiped corresponds to the center of a face-face contact between bipyramids 
of opposite orientations, accounting for all contacts between oppositely-oriented 
bipyramids. Four of these contacts are within one layer and four of them are 
with the layers above and below. The contacts formed between like-oriented 
bipyramids vary with the parameter x: for all values of x there are two edge- 
edge contacts, a vertex-edge contact and an edge-vertex contact (all of these 
contacts occur on two- fold axes and are within one layer); for x — 29/56 there 
are four additional edge-edge contacts with dimers in neighboring layers (which 
turn into overlaps for x < 29/56); and for x = 9/14 there are instead two 
vertex-face contacts and two face-vertex contacts with dimers in neighboring 
layers (which again turn into overlaps for x > 9/14). Thus, each dimer makes 
respectively twelve, sixteen, or sixteen contacts in the three cases, and corre- 
spondingly, each tetrahedron makes eight, ten, or eleven contacts. 



3 Simple double-lattice packing 

Our numerical search also yielded a packing with two tetrahedra per repeating 
unit, which we could identify as a simple double lattice (that is a double lattice 
of tetrahedra, not of dimers) with packing fraction 139 t j gg' ,/ ^ ~ 0.7194. We 
present therefore a second theorem. 

Theorem 2. There exists a packing of regular tetrahedra having packing fraction 
139 ~369 V/ ^" ' -^ l * s packing is periodic, with each unit cell of the lattice contain- 
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fundamental 


Tn = conv-fr,- 1 i = 1 2 3 41 


tetrahedron 


ri = [(433 - 86v / 10)a + (611 - 133\/T0)b+ 




(188 - 22\/T0)c]/246 




r 2 = [(111 - 30^T0)a + (93 - 75Vl0)b+ 




(-66-42VT0)cl/246 




r 3 = [(-85 - 28VlO)a + (13 - 29 v / 10)b+ 




(4 + 10Vl0)c]/246 




r 4 = [(179 - 106VT0)a + (427 - 101%/l0)b+ 




(-20 - 50ylT))c]/246 


other tetrahedron 


T\ = I(T ) = -T 


in the unit cell 




space group 


translations by a, b, c 


generators 


/ = inversion about 


packing fraction 


(139 + 40%/l0)/369 


coordinate basis for 


a= (l,-(13-4\/T0)/3,0) 


which tetrahedra are 


b = (-(4- v / 10)/3,3- VT6,-1) 


regular with unit edge 


c = (3- VlO, 1,(4- Vw)/3) 



Table 2: The construction of the simple-double-lattice packing in a general 
triclinic coordinate basis. 



Figure 3: Small portions of one layer and three stacked layers in the simple- 
double-lattice packing. 



G 



Name 



N Z Transitive 



Optimal lattice[8] 
Warp and weft [3 
Welsh[3] 

Simple double lattice 
Wagon wheels [4] 
Compressed wagon wheels [5] 
Disordered wagon wheels [6 J 
Quasicrystal approximant[7J 
Dimer double lattice 



0.7786 
0.7820 
0.8226 
0.8503 



18/49 w 0.3673 
2/3 w 0.6666 
17/24 « 0.7083 



100/117 g 0.8547 




1 

2 

34 
2 

18 

72 

314 

656 

4 



No 

8 to 11 Yes 



14 Yes 

10 Yes 

25.9 No 

19 Yes 

7.1 No 

7.6 No 

7.4 No 



Table 3: Some studied transitive and non-transitive packings of regular tetra- 
hedra with packing fraction 0, number of tetrahedra in the repeating unit N, 
and average number of contacts per tetrahedron Z where available. 

ing two tetrahedra. The group of isometries leaving the packing invariant is a 
crystallographic space group of type PI and acts transitively on the individual 
tetrahedra of the packing. 

The proof of this theorem proceeds equivalently to the proof of Theorem 
1. The vertices of the initial tetrahedron To are given in Table 2 in terms 
of a triclinic coordinate basis which is also given in Table 2. The group of 
isometries we use to construct the packing is generated by translation by the 
three vectors a, b, and c and by inversion about the origin. This is a space 
group of crystallographic type PI, with a point group of order 2 and a triclinic 
lattice [11] . Figure 3 shows a portion of the packing. 

The simple-double-lattice packing again has eight inversion centers per prim- 
itive cell at the vertices of a parallelepiped. However, in this case only five of the 
vertices are on the surface of the tetrahedron. Each tetrahedron in the packing 
is in contact with nineteen others. 

4 Discussion 

In Table 3, we compare the packings presented here to other studied packings 
of regular tetrahedra. Both packings are denser than the densest previously- 
reported transitive packing, a double lattice presented by Conway and Torquato 
(which we call the "warp- and- weft" packing due to the interweaving arrange- 
ment of its tetrahedra) [3], and the dimer double lattice is denser than any 
previously-reported packing. 

The results presented go against the recent trend of ever-growing repeating 
units in densest-known packings and demonstrate that a large repeating unit is 
not a necessary property of a dense packing of regular tetrahedra. It is curious 
that previous simulations, utilizing a more physical search dynamic [7J, 
yielded dense packings that were either disordered, had quasicrytalline order, or 
had crystalline order characterized by a very large repeating unit, and were not 
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able to find the denser class of structures presented here, (reminiscent perhaps 
of Kurt Vonnegut's ice-nine, a fictional phase of water that is more stable, but 
kinctically unreachable) . 

Our results yield the surprising situation wherein the densest-known packing 
of icosahedra is now sparser than the corresponding packing of tetrahedra, a 
solid which just four years ago was a prime candidate for a counterexample 
of a conjecture by Ulam that the sphere is the sparsest-packing convex solid 
[3]. As the packing can be generally extended to any tetrahedron in a three- 
parameter family generated by deformations of the monoclinic coordinate basis, 
if any tetrahedron provides a counterexample of Ulam's conjecture, it is not a 
tetrahedron of that family. 

The regular tetrahedron is no longer outcast, as it long was, from the re- 
spectable family of convex polyhedra whose largest-achieved packing density is 
realized by a transitive arrangement. While there are some convex solids whose 
maximum packing density clearly cannot be achieved by a transitive arrange- 
ment (the convex Schmitt-Conway-Danzer polyhedron can tile space, but only 
in aperiodic and non-transitive ways [5]), the majority of regular and semi- 
regular polyhedra have been to-date packed most densely in transitive packings 
[5j. Whether this situation is accidental, the result of bias favoring the discovery 
of transitive packings, or a more fundamental property governing the packing 
of a certain class of solids is still an open question. 

This work was supported by grant NSF-DMR-0426568. 
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